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ABSTRACT 

We present an atlas of SpitzerflRS high resolution (R^600) 10-37/im spectra for 24 well known starburst 
galaxies. The spectra are dominated by fine-structure lines, molecular hydrogen lines, and emission bands of 
polycyclic aromatic hydrocarbons. Six out of the eight objects with a known AGN component show emission 
of the high excitation [Ne V] line. This line is also seen in one other object (NGC 4194) with, a priori, no known 
AGN component. In addition to strong polycyclic aromatic hydrocarbon emission features in this wavelength 
range (n.3, 12.7, 16.4 fim), the spectra reveal other weak hydrocarbon features at 10.6, 13.5, 14.2 fim, and a 
previously unreported emission feature at 10.75/im. An unidentified absorption feature at 13.7/im is detected 
in many of the starbursts. We use the fine-structure lines to derive the abundance of neon and sulfur for 14 
objects where the HI 7-6 line is detected. We further use the molecular hydrogen lines to sample the properties 
of the warm molecular gas. Several basic diagrams characterizing the properties of the sample are also shown. 
We have combined the spectra of all the pure starburst objects to create a high S/N template, which is available 
to the community. 

Subject headings: galaxies: starburst. Infrared: galaxies, ISM: lines and bands 



1. EVITRODUCTION 

As galaxies evolve, they often go through periods when 
their star formation rates increase substantially. This is usu- 
ally ascribed to large-scale perturbations in their structure, 
such as the formation of bars triggered by interactions with 
companions, leading to the accumulation of large quantities 
of gas in their nuclei. A small number of those so-called 'star- 
burst' galaxies (e.g. Weedman et al. 1981) have been identi- 
fied in the local Universe and have been studied extensively at 
all wavelengths. At high redshifts, obscured starburst galax- 
ies are much more common (Blain et al. 2002; Elbaz & Ce- 
sarsky 2003), and are believed to give rise to the bulk of the 
cosmic infrared background radiation. Reviews of their prop- 
erties can be found in Kennicutt (1998); Lagache et al. (2005) 
and Lonsdale et al. (2006). 

Understanding the physics of starburst galaxies in the 
local Universe is of great importance, as they are the closest 
laboratories we have for studying the starburst phenomenon, 
and because they may represent 'scaled-down' versions of 
the starbursts seen at high redshift. The obscured nature of 
(most) starburst galaxies, regardless of their redshifts, how- 
ever means that probing their nuclear regions in the optical is 
challenging. This task is better done in the infrared, where 
dust opacity is lower. The advent of major infrared space 
facilities such as IRAS (Neugebauer et al. 1984), ISO (Kessler 
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et al. 1996) and Spitzer (Werner et al. 2004), and in particular 
the accessibility to high quality spectroscopy, has allowed 
great progress in this field, including detailed studies of the 
mid-infrared (MIR) region of the spectrum in galaxies. It was 
discovered that the mid-infrared (MIR) is rich in emission 
and absorption features of gas and dust. Emission from Poly- 
cyclic Aromatic Hydrocarbons (PAHs, Leger & Puget 1984) 
was shown to be a reliable star formation tracer (e.g. Roussel 
et al. 2001; Forster Schreiber et al. 2004; Peeters et al. 2004), 
though variations with metallicity have also been observed 
(Wu et al. 2006; Smith et al. 2007; Rosenberg et al. 2008). 
The use of infrared atomic lines with different degrees of ion- 
ization has proven invaluable in studying the properties of the 
ionization field in starbursts (e.g. Sturm et al. 2002; Verma et 
al. 2003; Wu et al. 2006; Dale et al. 2006). In ultra-luminous 
infrared galaxies (ULIRGs), diagnostic diagrams based on 
line ratios such as [NeVA14.32^m]/[NeIIA12.81^m], 
[OIV A25 . 89/im]/[NeII A 1 2.81 fim] , 

[SIIIA33.48/im]/[SiIIA34.82^m] and combinations of 
those with PAH emission features have been developed to 
reveal their nature (e.g. Genzel et al. 1998; Laurent et al. 
2000; Thornley et al. 2000; Dale et al. 2006; Armus et al. 
2007; Farrah et al. 2007, 2008). Additional diagnostics 
using the silicate absorption features (Spoon et al. 2007) or 
the strength of the rotational molecular hydrogen lines (e.g. 
Rigopoulou et al. 1999; Roussel et al. 2007) have also been 
proposed. 

As part of our continuing efforts to provide a comprehen- 
sive picture of the starburst phenomenon at low redshifts, 
we here present an atlas of high resolution (R'-^600) 9.7- 
37/im spectra of well known nearby starburst galaxies, taken 
with the Infrared Spectrograph^ (Houck et al. 2004) on-board 
Spitzer This follows on from our earlier work on the same 
sample using the low resolution IRS modules (Brandl et al. 
2006). We have created a high S/N template from these spec- 

' The IRS was a collaborative venture between Cornell University and Ball 
Aerospace Corporation funded by NASA through the Jet Propulsion Labora- 
tory and the Ames Research Center. 
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TABLE 1 
The Starburst Galaxy Sample 



Name 


Type 


RA 


DEC 


z 


D 


log[LiR] 


AOR ID" 


SH t;^ 


LH t,„, 












Mpc 


L© 




(cycles X s) 


(cycles X s) 


NGC253 


SB 


00''47'"33.12-' 


-25'' 17"" 17.6* 


0.0008 


2.5 


10.23 


9072640 


12 X 6 


12 X 6 


NGC520 


SB 


0l''24'"35.07'' 


+03''47""32.7* 


0.0076 


30.2 


10.91 


9073408 


4 X 30 


3 X 60 


NGC660 


SB+AGN'' 


0l''43'"02.35* 


+13''38""44.4* 


0.00283 


12.3 


10.49 


9070848 


4 X 30 


8 X 14 


NGC 1097 


SB+AGN'' 


02''46'" 19.08* 


-30^* 16""28.0* 


0.00456 


16.8 


10.71 


3758080 


4 X 30 


2 X 60 


NGC 1222 


SB 


03''08""56.74'' 


-02''57""18.5* 


0.00818 


32.3 


10.60 


9071872 


4 X 30 


2 X 60 


NGC 1365 


SB+AGN'' 


03''33'"36.37* 


-36''08""25.5* 


0.00545 


17.9 


11.00 


8767232 


4 X 30 


8 X 14 


IC342 


SB 


03'"46'"48.51' 


+68''05'"46.0* 


0.00010 


4.6 


10.17 


9072128 


4 X 30 


8 X 14 


NGC 1614 


SB 


04*33'"59.85' 


-08'' 34*" 44.0* 


0.016 


62.6 


11.60 


3757056 


4 X 30 


8 X 14 


NGC 2146 


SB 


06*18"'37.71' 


+78''21"'25.3* 


0.003 


16.5 


1 1 .07 


9074432 


1 X 6 


1 X 60 


NGC 2623 


SB+AGN'' 


08*38"'24.08' 


+25'' 45"' 16.9* 


0.01846 


77.4 


11.54 


9072896 


4 X 30 


2 X 60 


NGC 3079 


SB+AGN'' 


10*01"'57.80' 


+55'' 40*" 47.1* 


0.00372 


17.6 


10.52 


3755520 


4 X 30 


2 X 60 


NGC 3256 


SB 


10* 27"" 5 1.27' 


-43'' 54"! 13.8* 


0.00913 


35.4 


11.56 


9073920 


1 X 6 


1 X 60 


NGC 33 10 


SB 


10*38"'45.96' 


+53''30"'()5.3* 


0.0033 


19.8 


10.61 


9071616 


4 X 30 


2 X 60 


NGC 3556 


SB 


11* 11*" 30.97' 


+55'' 40*" 26.8* 


0.00233 


13.9 


10.37 


9070592 


4 X 30 


2 X 60 


NGC 3628 


SB+AGN'' 


11*20" 17.02* 


+13''35'"22.2* 


0.00281 


10.0 


10.25 


9070080 


4 X 30 


2 X 60 


NGC 4088 


SB 


12*05'"34.19* 


+50''32'"20.5* 


0.00252 


13.4 


10.25 


9070336 


4 X 30 


2 X 60 


NGC 4194 


SB+AGN" 


12*14'"09.64* 


+54''31"'34.6* 


0.00835 


40.3 


11.06 


3757824 


4 X 30 


2 X 60 


Mrk52 


SB 


12*25'"42.67* 


+00''34"'20.4* 


0.0071 


30.1 


10.14 


3753216 


4 X 30 


2 X 60 


NGC 4676 


SB 


12*46"" 10. lO* 


+30''43""55.0* 


0.022 


94.0 


10.88 


9073152 


8 X 30 


6 X 60 


NGC 48 18 


SB 


12*56'"48.90' 


-08''31""31.1* 


0.00355 


9.4 


09.75 


9071104 


4 X 30 


2 X 60 


NGC 4945 


SB+AGN'' 


13*05'"27.48* 


-49''28""05.6* 


0.00186 


3.9 


10.48 


8769280 


1 X 6 


1 X 60 


Mrk266 


SB+AGN'' 


13*38'" 17.69* 


+48''16"'33.9* 


0.02786 


115.8 


11.49 


3755264'' 


4 X 30 


2 X 60 


NGC 7252 


SB 


22*20™44.77* 


-24''40"41.8* 


0.0156 


66.4 


10.75 


9074688 


4 X 30 


3 X 60 


NGC 7714 


SB 


23*36"" 14. 10* 


+02''09""18.6* 


0.00933 


38.2 


10.72 


3756800 


4 X 30 


2 X 60 



'^This is the AOR ID of the high resolution observations only. 

''The evidence for an AGN in these objects come froin the literature (NED). 

"The evidence for an AGN comes from the detection of the [Ne V] line in this work. 

''The coordinates in the Table for this galaxy are slightly off the nucleus (see §2 and Fig. 1). The galaxy is still included in the LH slit but for the smaller 
SH slit we used an observation from a different program (PID=3237, AORkey=10510592) where the source is well centered. 



tra. Our goal is to make this set of high quality spectra avail- 
able to the community, and to report on their characteristic 
features. Observations are described in §2, and data reduction 
and analysis is described in §3. We present results in §4 and 
discussion in §5. Finally, we summarize our conclusions in 
§6. 

2. OBSERVATIONS 

Our sample is presented in Table 1 where information on 
the galaxy optical classifications, distances, magnitudes and 
other basic information on the observations is given. In nine 
objects there is evidence of an active galactic nucleus (AGN), 
even though it does not dominate the bolometric luminosity 
of the system. The sample was observed as part of the IRS 
Guaranteed Time program (PID 14) and consists of 24 well 
known nearby starburst galaxies, with an average distance of 
33 Mpc and an average infrared (IR) (8-1000/im) luminosity 
of 5x10'" L0. The observations were made using both the 
high- and low-resolution modules of the IRS. Here we present 
the high-resolution spectra. The low-resolution spectra can be 
found in the companion paper by Brandl et al. (2006). 

The two high-resolution modules, short-high (SH) and 
long-high( LH), cover the wavelength region between 10 
to 37 /im at a resolution of 600 and have slit sizes of 
4.7" X 1 1.3" and 1 1. 1" x 22.3" respectively. The data were 
taken using the Staring Mode observing template which pro- 
duces spectra in two nod positions. The separation between 
the nod positions is 3.8" for SH and 7.4" for LH. The ob- 
servations were centered on the nucleus of the galaxies (see 



Fig. 1) at the coordinates given in Table 1. However, as can 
be seen in Figure I, the observations of NGC 520 and Mrk52 
are slightiy mispointed with the SH slit measuring part of the 
nucleus. In NGC 3310 the mispointing is more severe and 
affects both modules and therefore any information we give 
about this galaxy in the paper refers to the position we mea- 
sured and not the nucleus. Mrk266 is also mispointed but 
the LH module is large enough to contain the whole galaxy 
and for the SH module we made use of an observation of an- 
other Spitzer program in which the source is well centered 
(PID=3237, AORkey=10510592). The area covered by the 
SH slit is equivalent to a spatial scale of about 0.085 kpc to 
3.2 kpc, depending on the distance to the objects in our sam- 
ple. No dedicated background observations were made. This 
means that there will be a small contribution of the zodia- 
cal light to the spectra in Figures 2, 3 and 4. This has no 
consequence for the line measurements. Foreground emis- 
sion of PAHs from our MW is also possible but this is much 
fainter than the targets themselves. For information of the 
mid-infrared continuum properties the reader is referred to 
the paper on the low resolution spectra (Brandl et al. 2006) 
where the background contribution to the continuum has been 
subtracted from the spectra. 

3. DATA REDUCTION & ANALYSIS 

The data were processed using version sl5.3 of the Spitzer 
Science Center's pipeline which is maintained at Cornell and 
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Fig. 1 . — IRAC 8 fim, and 2MASS-K images of the targets with the SH and LH slits overlayed on them. The grey scale of the image is linear up to the 50% 
level of brightness and is black above it. 



using a script version of Smart (Higdon et al. 2004). The 
reduction begins from the droop images which are equivalent 
to the more commonly used bed data but without flatfield and 
stray-cross-light corrections. Unstable sensitive pixels (usu- 
ally referred to as 'rogue' pixels) were replaced using the SSC 
tool irsclean^^, which uses campaign-based pixel masks to 
identify the bad pixels. At the location of a pixel identified as 
bad, the empirical profile in the cross dispersion direction is 
determined using the average profile in the neighboring rows. 
This profile is then scaled to the profile of the row containing 

Smart can be downloaded from this web site 
http://ssc.spitzer.caltech.edu/archanaly/contributed/smart/index.html. 

" This tool is available from the SSC web site: 
http://ssc.spitzer.caltech.edu 



the bad pixel using the remaining good pixels of that row. The 
rogue pixels are then replaced by the value given by the empir- 
ical profile. This method of replacing bad pixels is equivalent 
to interpolation when the bad pixel is located on the wings of 
the spatial profile but it has a definite advantage over interpo- 
lation when the bad pixel is located at the maximum of the 
spatial profile. 

For each module the individual observations were aver- 
aged together, and spectra extracted using full aperture extrac- 
tion. Flux calibration was performed by dividing the averaged 
spectra by the spectrum of the calibration star ^ Dra (extracted 
in the same way as the target) and multiplying by its theoreti- 
cal template (Cohen et al. 2003, Sloan et al. in prep). Fringes, 
when present, were removed using the IRS de-fringing tool 
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Fig. 2. — SH and LH spectra. The main lines and features are indicated in tlie top panel. The discontinuity around 19.5 /xm is caused by the extended nature 
of the sources combined with the different sizes of the SH and LH apertures, see also Brandl et al. (2006), their Figure 1. 
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Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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TABLE 2 

SH LINE fluxes". 



Object [SIV] H2(S2) HI(7-6) [Nell] [Ne V] [CI II] [Ne III] HjCSl) [PIII] [Fell] [S III] 
10.51'' 12.28 12.37 12.81 14.32 14.38 15.55 17.03 17.89 17.95 18.71 



NGC 253 


<10.50 


64.32 


<27.45 


2832.33 


<20.5C 






6.78 




64.20 




NGC 520 


<0.47 


3.49 


<0.65 


44.62 


<0.64 






0.22 




1.25 




NGC 660 


2.79 


12.53 


<3.68 


353.01 


<2.26 




0.47 


1.01 




7.55 




NGC 1097 


<0.63 


6.38 


0.57 


37.24 


<0.43 






0.50 


0.11 


9.27 




NGC 1222 


22.24 


2.93 


2.12 


80.57 


<0.57 




0.45 


0.16 


0.20 


1.30 




NGC 1365 


23.37 


9.00 


3.18 


139.50 


18.83 




1.42 


0.91 


0.39 


3.27 


0.51 


IC342 


4.76 


10.32 


6.22 


615.46 


<2.41 




0.66 


0.96 


0.79 


10.52 




NGC 1614 


6.89 


5.14 


1.77 


249.00 


<0.99 




0.54 


0.51 


0.37 


7.00 




NGC 2146 


6.30 


12.79 


6.81 


625.00 


<2.81 




0.44 


1.35 


1.53 


15.35 




NGC 2623 


1.20 


3.15 


<0.49 


55.55 


2.71 




0.15 


0.19 




0.66 


0.15 


NGC 3079 


<0.96 


23.81 


<2.09 


98.96 


1.31 






1.28 




2.03 


0.21 


NGC 3256 


5.25 


15.04 


6.19 


514.19 


<2.35 




0.57 


1.07 


1.38 


9.34 




NGC 3310 


4.46 


<0.56 


0.58 


27.57 


<0.25 




0.41 




0.13 


0.61 




NGC 3556 


<0.56 


1.06 


0.71 


21.47 


<0.37 






0.16 


0.18 


0.42 




NGC 3628 


<0.90 


6.53 


<1.68 


125.63 


0.90 






0.43 




2.37 


0.35 


NGC 4088 


<0.67 


2.70 


0.95 


37.03 


<0.39 






0.32 


0.11 


1.16 




NGC 4194 


11.46 


5.03 


2.12 


165.47 


2.96 




0.75 


0.51 


0.58 


3.90 


0.65 


JVLTK JZ 


i.uz 


i.jy) 


<U.ol 


OQ IQ 


<U.D / 




0.10 


0.22 




0.53 




NGC 4676 


0.56 


2.03 


0.48 


27.78 


<0.25 




0.06 


0.14 


0.09 


0.96 




NGC 48 18 


2.21 


5.87 


1.71 


184.93 


<1.21 




0.16 


0.35 


0.50 


5.60 




NGC 4945 


0.72 


33.93 


<3.66 


583.84 


3.84 




0.18 


0.85 




12.50 


0.34 


Mrk266 


9.00 


3.72 


0.69 


57.04 


7.96 




0.38 


0.32 


0.06 


1.64 


0.17 


NGC 7252 


<0.56 


2.70 


<0.57 


41.68 


<0.33 






0.29 




1.94 




NGC 7714 


14.76 


2.85 


2.18 


102.55 


<1.00 




0.63 


0.41 


0.36 


1.98 




Template"^ 


14.16 


13.86 


5.56 


452.83 





25.57 


204.64 


67.24 


22.89 


25.24 


666.37 


5.92 


9.59 


2.35 


6.01 


5.97 


14.93 


0.89 


7.53 


10.82 


<0.51 


<0.48 


9.10 


0.15 


0.21 


0.35 






0.59 


4.05 


36.96 


23.52 


4.89 


2.52 


90.02 


0.72 


0.49 


0.90 


1.00 


0.36 


1.71 


<0.47 


6.33 


13.22 


<0.99 


<0.51 


15.88 




0.18 


0.36 






0.40 


<0.88 


89.41 


9.53 


1.49 


0.83 


64.80 




1.74 


0.31 


0.21 


0.26 


1.13 


2.05 


59.53 


18.42 


2.22 


1.72 


53.72 


0.37 


0.92 


0.63 


0.38 


0.38 


3.58 


7.99 


37.20 


11.69 


8.27 


6.55 


320.03 


0.77 


0.90 


0.76 


0.79 


0.84 


5.96 


1.15 


63.32 


9.42 


5.14 


2.13 


83.03 


0.32 


1.59 


0.58 


0.58 


0.51 


2.64 


5.48 


91.16 


26.12 


7.71 


3.24 


190.12 


0.66 


0.95 


1.64 


0.42 


0.63 


4.15 


0.69 


15.08 


7.27 


<0.52 


<0.52 


8.20 


0.13 


0.24 


0.23 






0.36 


1.57 


23.17 


40.63 


<0.52 


0.90 


11.83 


0.21 


0.46 


1.04 




0.13 


0.16 


5.99 


64.42 


27.99 


6.29 


3.69 


171.83 


0.85 


0.86 


2.84 


0.73 


0.53 


2.07 


0.29 


28.35 


1.88 


<0.61 


<0.59 


20.50 


0.10 


0.45 


0.11 






0.31 


<0.40 


3.23 


2.90 


<0.45 


<0.34 


12.36 




0.09 


0.12 






0.22 


2.35 


10.05 


16.24 


<0.79 


<0.96 


21.32 


0.18 


0.29 


0.60 






0.77 


0.51 


2.45 


4.72 


<0.40 


<0.37 


12.50 


0.16 


0.14 


0.16 






0.26 


2.49 


53.99 


9.25 


1.70 


1.90 


70.92 


0.38 


0.78 


0.66 


0.32 


0.38 


1.44 


<0.43 


3.82 


1.80 


<0.45 


<0.39 


28.45 




0.29 


0.19 






1.63 


<0.31 


4.68 


5.19 


<0.36 


<0.36 


11.32 




0.16 


0.32 






0.37 


2.27 


13.73 


16.83 


2.61 


1.41 


71.96 


0.35 


0.45 


0.57 


0.35 


0.18 


1.05 


6.51 


68.98 


54.15 


1.77 


2.85 


30.53 


0.27 


14.25 


4.22 


0.25 


0.24 


0.58 


0.46 


27.95 


8.40 


0.67 


0.48 


24.33 


0.07 


0.76 


0.97 


0.15 


0.14 


0.69 


0.56 


3.72 


5.10 


<0.54 


<0.57 


11.96 


0.10 


0.14 


0.18 






0.31 


<0.68 


77.42 


4.43 


2.01 


1.58 


81.50 




0.92 


0.35 


0.25 


0.27 


1.10 


4.58 


100.00 


30.83 


5.74 


3.30 


187.97 



^Per object, the line fluxes are given in the first row in units of 10"^^ W cm"^, with the uncertainties shown in the row below. Note that, except 
for the values of the template, the fluxes listed here are not to be directly compared (unless scaled) to those in Table 3 because of the different 
aperture used. 

^ Rest wavelength in ^m. 
The line fluxes for the template spectrum are relative to F([Ne III])=100. 



irsfringe^^. Finally, the spectra for each nod position were 
combined. It is worth mentioning that in six of the star- 
bursts (IC 342, NGC 520, NGC 1097, NGC 1614, NGC 3310, 
and NGC 3628) the continua differed between nod positions. 
These sources are extended (see Figure 1), and this difference 
may reflect the different regions the nod positions are sam- 
pling. Since we are interested in the overall properties of the 
emission lines and PAH bands, and differences in these quan- 
tities lie within the errors in each nod position, we combined 

This contributed software can be downloaded from 
http://ssc.spitzer.caltech.edu/archanaly/contributed/irsfringe/ 



the nod positions in these starbursts as it was done for the rest 
of the sample. 

Fine- structure Une fluxes are derived by fitting a Gaussian 
profile and integrating the flux above a local continuum. A 
polynomial of order one was used to fit the local continuum 
for a given line, except for the [Ne II] line which sits at the 
sharp flank of a PAH feature (see Fig. 3). For this line, a third 
order polynomial was used to fit the PAH underneath the Une. 
This is done for each nod separately and then the flux in both 
measurements is averaged. Errors bars were estimated using 
the difference between the two nod positions or the standard 
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TABLE 3 

LH LINE FLUXES^. 



Object 


rc^ TTTI 

[re lilj 
22.93'' 


rxT^ \7'\ 

[JNe VJ 
24.31 


[OIVJ 
25.89 


rt:^ TTI 

25.98 


H2(S0) 
28.22 


rc TTTI 

[S lilj 
33.48 


rc; TTI 

[Si ilj 
34.81 


rxT^ TTTI 

[Ne 1111 
36.01 




120.84 


<73.36 


154.74 


250.56 


<77.55 


1 538 03 


2412.03 


-243.86 




18 17 




26.95 


24 43 




30.06 


48 02 






5 95 


<1 42 


8 10 


7.56 


8.23 


89.44 


190 71 


<6 90 




0.67 




0.88 


78 


0.60 


8.60 


3 41 






12 33 


<6 59 


18 80 


22 92 


10 58 


246 07 


441 96 


<8 96 




1 64 




2 57 


1 82 


2 67 


3 21 


3 56 




IN v) V. 1 yjy 1 


4 20 


<1 07 


4 41 


1 1 48 


6 59 


100 05 


225 19 


<5.29 




0.18 




1 21 


1 07 


69 


3 95 


7 00 






7.04 


<0 78 


9 92 


4 74 


3 25 


132 52 


1 12 05 


12 72 




0.38 




48 


38 


34 


2 30 


1 70 


1 30 








141 57 


22 20 


15 91 


246 89 


500 30 


< 14 78 




1 20 


1 18 


2 06 


2 33 


1 51 


12 17 


17 78 




IC 342 


36 43 


<4.95 


<7 70 


51 87 


<9 10 


672 46 


985 73 


<29 96 




2.69 






1.89 




7.28 


10.29 




iWJK^ ID It 


13 14 


97 


8 A8 


19 QQ 




ini ofi 


If O.DW 


12 15 




2 10 




U.oj 


1 1 S 
1. 1 J 




2 11 


4 11 


3 44 






<^'^ 70 


17. J.J 


S9 ftO 


ID.DD 


09 








2.58 




4.48 


4.59 


2.33 


36.93 


21.70 




Nnr 969'^ 


<1.31 


2.21 


9.62 


2.23 


<2.08 


13.78 


28.74 


<9.41 






0.14 


0.72 


0.51 




1.41 


1.89 




NOP "^070 








14 57 


1 1 02 


56 88 


173 51 


^7 OS 




0.21 




0.84 


0.62 


0.75 


2 70 


5 83 






31.55 


<3.85 


12 23 


33.56 


12 70 


484 64 


623 37 


<33.15 




2 77 




3 02 


1 S9 


1 '^9 


1 '\ 7Q 


9 81 






5 66 


<1.00 


4 03 


8.78 


<2.32 


106 02 


143.54 


8 01 




0.59 




0.64 


0.74 




1.17 


1.58 


1.19 


Nnr '^ssfi 


2 24 


<1.18 


<1.50 


3.14 


4.55 


96 37 


96.44 


<3.00 




0.34 






47 


47 


1 18 


80 




Nnr '^fi9s 




^1 0^ 




1 1 96 


1 1 57 


149 93 


277 55 


1 OQ 




0.60 






1 07 


42 


2 20 


7.75 




Nnr 40X8 


0.98 


<r(\ SO 


0.73 


2.15 


3.38 


35.02 


63.33 


<1.33 




0.09 




0.21 


0.27 


0.23 


0.41 


0.73 




Nnr 41Q4 


9.96 


4.04 


27.45 


8.48 


<2.66 


151.01 


185.82 


<10.83 




1.00 


0.23 


1.63 


0.96 




1.65 


3.31 




Mrk52 






1 '^1 
i.ji 


1 on 


<^ 1. 




'^R 09 


<2 44 




0.61 




0.23 


0.31 




1.07 


0.60 




Nnr 4fi7fi 


<0.58 


<0.33 


1 37 


1 78 


1.89 


23.95 


41.27 


<2.49 








0.23 


0.22 


0.25 


1.51 


1.42 




NGC 4818 


7.96 


<1.96 


<1.59 


10.83 


4.68 


80.54 


128.66 


<5.61 




0.50 






0.88 


0.51 


2.44 


1.97 




NGC 4945 


8.69 


<7.95 


47.85 


55.21 


44.56 


338.91 


846.23 


-160.25 




2.09 




6.95 


5.35 


8.72 


70.93 


50.98 




Mrk266 


2.15 


11.12 


52.94 


4.35 


<1.60 


50.93 


87.02 






0.29 


0.59 


5.66 


0.70 




2.91 


4.44 




NGC 7252 


1.36 


<0.51 


1.33 


3.12 


1.10 


24.75 


70.53 


<2.20 




0.26 




0.25 


0.20 


0.33 


0.47 


0.43 




NGC 7714 


7.88 


<1.77 


5.51 


5.93 


<1.84 


115.90 


102.97 


7.45 




0.70 




0.94 


0.54 




2.45 


1.63 


1.29 


Template"^ 


11.25 




7.85 


16.07 


9.51 


267.16 


364.23 


8.10 



^For each object, the line fluxes are given in the first row in units of 10""^' W cm"^, with the uncertainties given 
in the row below. Note that the fluxes hsted here are not to be directly compared (unless scaled) to those in Table 2 
because of the different aperture used. 

''Rest wavelength in iim. 

*^ The line fluxes for the template spectrum are relative to F([Ne ni])=100. 



deviation of the residuals of the fit, whichever is the larger. 
Upper-limits are derived by measuring the flux of a Gaussian 
with a height three times the local rms, and with a FWHM 
equal to the instrumental resolution. Note that, except for the 
line measurements of the template. Tables 2 and 3 quote the 
line fluxes without any normalization applied to correct for 
the difference in aperture between the SH and LH modules. 

Following Hony et al. (2001); Vermeij & van der Hulst 
(2002);and Peeters et al. (2002), we have measured the 
strength of the PAH emission features by integrating the flux 
of the PAH bands above an adopted spline interpolated local 



continuum. For a given band the anchor points used to de- 
fine the continuum (see Table 4) are the same for all sources, 
as are the integration limits. Spectral Unes superimposed on 
the PAH bands were removed before integration. Our spline 
method for measuring PAH fluxes differs from that of PAH- 
FIT (Smith et al. 2007) and others (e.g. Boulanger et al. 1998; 
Li & Draine 2001) by adopting a continuum (which ignores 
the plateau) and integrating under the bands, rather than as- 
suming a PAH band profile and performing a simultaneous fit 
to all the bands. Galliano et al. (2008) concluded that both 
methods reveal similar trends in the data. The uncertainties in 
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TABLE 4 

PAH FLUXES^''. 



Object 


F(10.6)'' 


F(10.7f 


F(ll.O) 


F(11.3) 


F(12.0) 


F( 12.74) 


F(13.5) 


F(14.2) 


F( 16.45) 


F( 16.74) 


F(17.0) 






3.39 


19.00 


405.27 


18.67 


257.50 


19.85 


13.78 


23.67 




15.86 






0.35 


1 .24 


20.26 


2.30 


12.88 


2.16 


3.30 


4.98 




6.98 




23 


08 


91 


17 21 


1 1 1 


9 66 


34 


32 


86 


03 


1.36 




02 


02 


05 


86 


06 


48 


06 


02 


13 


0.01 


0.07 


i'fVJw uuu 


0.52 


77 


6.07 


83.40 


4 60 


43.45 


1 05 


2 09 


5.64 




4.06 




21 


0.10 


30 


4 17 


44 


2 17 


26 


22 


32 




0.38 




30 




1 17 


18 84 


99 


6 32 


29 


32 


1 06 




0.99 




02 




06 


94 


05 


32 


07 


02 


23 




0.05 




0.66 


0.26 


2 72 


27 56 


1 09 


9 69 


38 


27 


1 54 


33 


1.53 




0.12 


0.08 


14 


1 38 


18 


48 


10 


09 


17 


13 


22 




70 


0.51 


3 64 


51 38 


2 45 


20 04 


1 14 


65 


2 64 


31 


3 20 






98 


46 


9 S7 
Z.J / 


49 


1 08 


23 


27 


32 


23 


46 


IC342 


1.41 


0.64 


9 70 


1 19 19 


6 38 


47 77 


3.04 


3 26 


6 21 




4 67 




0.16 


0.1 1 


49 


5.96 


0.32 


2.39 


0.15 


16 


45 




23 




1.18 


0.12 


3.82 


47.29 


2.46 


22.80 


0.30 


1.98 


3.35 


0.34 


3.93 




u.uo 


U.Uj 






12 


1 14 


n 1 1 


n in 


60 


A A'^ 


n 9A 




1. JO 


1. 


IM-.OU 


168 24 


7 0*^ 


Rd. S7 


1. Jo 






1 10 

1.17 


12 16 




0.21 


0.13 


0.74 


8.41 


0.41 


4.24 


0.31 


0.36 


0.50 


0.36 


0.75 




0.29 


0.15 


0.91 


12.61 


0.63 


7.37 


0.43 


0.40 


1.13 




0.74 




0.04 


0.02 


0.05 


0.63 


0.06 


0.37 


0.03 


0.03 


0.06 




0.07 








% 99 


56 68 


3 13 


35 39 


1 58 


1 59 


5 07 


24 


3 41 








U. it) 


2 83 


16 


1 77 


08 


08 


25 


02 


17 




1.31 


0.86 


8.65 


108 19 


5 41 


50 18 


1 50 


1 69 


7 77 


99 


6 27 




0.31 


0.17 


43 


5 41 


30 


2 51 


12 


13 


43 


29 


0.55 








1 21 


14 79 


48 


5 03 


44 


31 


78 


A AR 


A 74 








06 


74 


02 


25 


02 


02 


10 


0.01 


0.04 








84 


8.46 


0.46 


3 17 


31 


29 


0.44 




0.44 








0.04 


42 


02 


16 


05 


01 


Oil 




0.03 




02 


23 


1 81 


29 98 


1 48 


19 88 


01 


63 


2 04 


19 


2 38 




02 


01 


09 


1 50 


07 


QQ 


14 


VJ.VJJ 


A 97 


A A1 


A 1 9 








1 21 


12 64 


w.DD 


S flQ 

■j.yjy 


n 17 


22 


\J.yj 




A ^iR 








0.06 


0.63 


0.03 


0.25 


0.01 


0.02 


0.09 




0.03 




0.88 


0.27 


4.84 


49.18 


1.55 


20.15 




1.00 


2.83 


0.31 


3.19 




0.14 


0.09 


0.24 


2.46 


0.21 


1.01 




0.16 


0.19 


0.15 


0.27 


Mrk52 








J.oo 




9 4.Q 




n 17 


A 

U.J J 




n '^1 

U. Jl 










90 


n ns? 

\j.\jO 


12 




U.UJ 


A A4 




A A'^ 

U.UJ 




0.23 




0.79 


9.41 


0.35 


4.18 


0.09 


0.43 


0.68 


0.05 


0.68 




0.02 




0.04 


0.47 


0.03 


0.21 


0.02 


0.02 


0.03 


0.02 


0.06 


NGC4818 


0.59 


0.48 


3.72 


43.53 


2.07 


18.95 


0.55 


1.18 


3.14 


0.41 


1.87 




0.11 


0.03 


0.19 


2.18 


0.12 


0.95 


0.09 


0.06 


0.16 


0.10 


0.17 


NGC4945 






1.36 


54.64 


3.64 


83.84 


4.85 


2.64 


6.20 


0.17 


3.49 








0.13 


2.73 


0.57 


4.19 


0.24 


0.13 


0.31 


0.06 


0.17 


Mrk266 






0.27 


3.71 


0.27 


1.07 




0.19 


0.23 




0.34 








0.06 


0.19 


0.11 


0.22 




0.06 


0.06 




0.13 


NGC7252 


0.14 




2.05 


23.54 


1.12 


7.55 




0.50 


1.31 




1.54 




0.04 




0.10 


1.18 


0.06 


0.38 




0.05 


0.07 




0.08 


NGC7714 


0.52 


0.27 


3.18 


28.21 


1.53 


11.34 


0.22 


0.50 


1.83 


0.21 


1.81 




0.14 


0.05 


0.16 


1.41 


0.15 


0.57 


0.07 


0.09 


0.11 


0.09 


0.16 


Template^ 


1.24 


0.59 


8.57 


100.00 


4.65 


41.26 


0.61 


1.31 


6.58 


0.47 


5.04 



"PAH fluxes and uncertainties (in the row below the measurement) in units of 10 ™W cm ^. 

''The anchor points for the continuum were selected at these positions: 10.21, 10.38, 10.64, 10.83, 11.74, 12.21, 13.08, 1 3. 47 ;im for PAHs between 10 
and 13Atm; 13.01, 13.29, 14.55, 14.70/im for PAHs between 13-15/jm; and 15.64, 15.83, 16.10, 16.61, 16.90, 17.70/jm for PAHs between 16-18/im. The 
integration limits for the PAHs (in the same order as they are in the table) are: 10.42-10.66, 10.66-10.83, 10.85-11.1, 11.1-11.7, 11.74-12.20, 12.20-13.1, 
13.4-13.68, 13.9-14.1, 16.2-16.62, 16.62-16.88, and 16.9,17.55(um. 

'^PAH position in parenthesis in /jm. 

''The identification of this feature is not know, but for simplicity we fold its measurement into this table. 
^The fluxes of PAH features in the template have been normalized so that F(l 1.3)=100. 



the PAH fluxes (Table 4) are based on the differences in the 
fluxes measured from the two individual nod spectra. How- 
ever, this does not take into account all the systematic uncer- 
tainties (especially flux calibration). We estimate these to be 
5% and choose it as our minimum uncertainty. 

4. RESULTS 

4.1. Individual Spectra 

The spectra are shown in Figure 2. The files used to gen- 
erate this plot are available in electronic format (tabl-24.txt 



files). Measurements of fine-structure lines and PAH features 
that fall within the wavelength range of the IRS SH and LH 
modules are given in Tables 2, 3, and 4 respectively. 

The spectra contain a variety of emission features. All 
of the spectra contain the prominent fine-structure lines of 
[Nell] (12.81/im), [Ne III] (IS-SS^m), and [S III] (18.71 and 
33.48/im). Higher ionization lines such as [O IV] (25.98/im) 
and [S IV] (10.5 1/xm) are seen in most of the sample, but 
are weak. The [Ne V] line at 14.32/im is detected in seven 
sources. The [Ar V] line at 13.1/um is only detected in 
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Fig. 3. — Template spectrum created by combining the Starburst galaxies in our sample. The grey area is the standard deviation in the mean indicating the 
spectral diversity of the starburst sample. 



Mrk266. In addition, many objects show emission of [Fe II] 
(25.99/im) and [Fe III] (22.92/xm). The fine-structure lines in 
the IRS spectra are unresolved (even those with AGN), except 
for Mik 266 and NGC 3079 where there is evidence of broad- 
ening in the [Ne II] and [Ne ni] lines (Spoon et al. in prep.). 
The recombination line HI (6-7) at 12.37/im is detected in 14 
objects. Rotational lines of H2, with the strongest transition 
0-0 S(l) emitting at 17.03^m, are present in most of the spec- 
tra. A number of PAH features are seen, the most prominent 
being the neutral C-H out-of-plane bending bands at 11.2/xm 
and the 12.7/im. The 16.4/7,m feature and 17/im complex are 
also detected throughout the sample. The high S/N allow us 
to detect in the spectra other weaker PAH features and a few 
'unidentified' bands. These are best seen in Figures 3 and 4, 
and are discussed in the next section. 

4.2. Starburst Template 

We have constructed a high-resolution starburst template 
spectrum from the 15 spectra in our sample with no signs 
of AGN activity (see Table 1). The template is available in 
electronic format (tab25.txt). This is equivalent to the star- 
burst template derived by Brandl et al. (2006) but using the 
IRS high-resolution modules. To create this template, the SH 
spectrum of each object was first scaled up to the LH spectrum 
using the overlapping continuum around \9iim. Based on a 
spline-defined continuum (with anchor points given in caption 



of Figure 4) each spectrum were separated into a continuum- 
subtracted spectrum (feature spectrum) and the associated 
continuum spectrum (featureless spectrum). The continuum 
template was generated by averaging every individual contin- 
uum spectrum after normalization to the 14/Ltm flux. Like- 
wise, the feature template spectrum was created by comput- 
ing the noise-weighted mean of every individual continuum- 
subtracted spectrum, after normalization to the total flux of 
the emission features. Finally, the starburst template spec- 
trum was assembled by scaling the continuum spectrum to 
the average ratio between the sum of the PAHs and line fluxes 
over the 14/im continuum before adding it to the feature tem- 
plate spectrum. The template is shown in Figure 3. Up until 
now, the 2-40/Ltm/5O-SWS spectrum of M&2 has been usually 
adopted in the Uterature as a starburst high-resolution tem- 
plate. However, M82 is a dwarf starburst with a low luminos- 
ity (3x10^" L0) while our template is built from a sample of 
15 galaxies spanning a range in IR color and luminosity typ- 
ical of starburst galaxies. In fact, seven of them are LIRGS 
(>10^' Lq), which in terms of galaxy evolution are much 
more important as they play a major role in the total energy 
production of the universe at z^l (Le Floc'h et al. 2005). The 
template we created may hence be more representative of the 
average physical parameters of local starburst galaxies, and 
the average line ratios and relative PAH strengths may serve 
as comparison data for studies of starburst galaxies at higher 
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Fig. 4. — Close-ups of the continuum subtracted template spectrum. The anchor points used to remove the continuum in the individual spectra to create this 
template are those used for the PAH integration (see caption of Table 4) plus the following points: 18.4, 19.4, 21, 22.5, 23.2, 25.0, 26.5, 28.0, 30.0, 32.0, 35.5, 
and 36.5/Lim 



redshifts. 

Figure 4 shows a zoom into several interesting parts of 
the template spectrum. Panel a shows that, in addition to 
the [S IV] line at 10.5 l/xm, the spectrum contains the rarely 
detected 10.60/im PAH feature and another broad feature at 
10.74/xm. Panel b shows details of the 11.3/xm PAH feature 
including the 11.0/im component. Sloan et al. (1999) noted 
that this feature (11.0/im) indicates ionized PAHs based on 
spatially resolved spectra of the reflection nebula NGC 1333 
SVS 3. Another PAH feature at 12.0/Ltm is seen followed by 
the H2 0-0S(2) line and the HI Hua recombination line at 
12.37/im. The [Ne II] line (12.8 1/im) sits in the flank of 
the 12.7/im PAH feature. In panel c there is a PAH feature 
at 13.53/im followed by an unidentified absorption feature 
(13.65-14.0/im) which has, to our knowledge, not been previ- 
ously detected. This absorption is detected in most of the star- 
bursts but the strongest absorptions are found in those sources 
also harboring an AGN. This is shown in Figure 5, where the 
objects have been ordered according to the strength of the ab- 
sorption feature, and where out of the bottom seven objects 
(those with stronger absorption), six have an AGN. Its detec- 
tion in NGC 4945 is also confirmed from a map study in the 
region where the dip is seen to change in strength at differ- 
ent positions in the central region of the galaxy (Spoon et al., 
in prep). While not reported, this dip is also present in the 



low-resolution starburst template (made from the same sam- 
ple) in Brandl et al. (2006). The weak 14.22/im PAH feature is 
present and the line at 14.36/jm corresponds to [CI II]. Panel d 
displays more PAH features at 16.4/im, 16.7/im, and 17.0/um. 
The H2 0-0 S(l) transition at 17.0 fj,m is strong and the [Fe II] 
line at 17.98/xm, while weak, is still detected. Another iron 
line ([Fe III] at 22.93/im) is shown in panel e. The last panel 
shows a blend of lines which correspond to [O IV] (25.89/;im) 
and [Fe II] (25.98/im), the stronger of which is the iron line. 

5. DISCUSSION 

5.1. Conditions in the Ionized Gas 

The presence of individual lines can set constraints on the 
ionization conditions in the hot gas in our sample. The most 
notable of these is [Ne V]. The [Ne V] has a high ionization po- 
tential (IP) of 97.1 eV and has never been detected in purely 
stellar environments, but is routinely detected in spectra of 
galaxies hosting an AGN (Sturm et al. 2002; Lutz et al. 2003; 
Weedman et al. 2005; Dale et al. 2006) and in about half of 
the 53 ULlRGs analyzed by Farrah et al. (2007) using Spitzer 
data. This line is also seen in young SNe remnants (Oliva 
et al. 1999) but in starbursts, where the overall shape of the 
mid-IR continuum is not flat, the presence of [NeV] would 
imply a weak AGN rather than SNe emission (the lifetime of 
the latter is too small to allow for a substantial contribution to 
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Fig. 5. — Continuum divided spectrum to liigliliglit (in grey) tlie absorp- 
tion feature between 13.56 and 14.0^m. Tlie sources liave been ordered by 
decreasing strength of the feature. The vertical lines indicate the position of 
the [Ne V] and [CI II] fine-structure hues. 



the spectram). From the objects Usted in Table 2 previously 
known to have an AGN component in the literature, the [Ne V] 
line is present in Mrk266, NGC 1365, NGC 2623, NGC 3079, 
NGC3628, and NGC 4945 (see also Dudik et al. 2007), but 
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Fig. 6. — Hardness of the radiation field as traced by the sulfur 
(10.51/ini/18.71/im) and neon (15.55/ini/12.81/im) line ratios. 
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Fig. 7. — Hardness of the radiation field as given by the [Ne III]/[Ne II] 
and [S IV]/[Ne II] ratios. The dashed line indicates the best fit to the data 
excluding NGC 4945: log([Afe///]/[iVe//]) = -0.38(±0.12)+0.68(±0.08) x 
\og([SIV]l[NeII]). 



not in NGC 660 and NGC 1097. In these two galaxies the 
[Ne V]14.3/xm/[Ne II]12.8/im hne ratio is <0.07 hinting at an 
AGN contribution of less than a few percent using the diagram 
by Farrah et al. (2007) (their Fig. 16). NGC 1097 is a partic- 
ularly interesting case; it also does not show [O rv] or [S IV], 
hence the AGN is likely very weak in the mid-lR. It either has 
a complex extinction pattern around it or it is overwhelmed by 
the starburst. We further detect [Ne V] in NGC 4194, which 
has no previous AGN classification. 

The strongest detections of [Ne V] relative to [Ne n] in 
our sample are found in Mrk266 and NGC 1365, at ratios of 
0.140±G.005 and 0.135±0.005, roughly a factor 10 above the 
rest of the sample. Adopting the AGN-starburst mixing sce- 
narios of Sturm et al. (2002) and Farrah et al. (2007), this 
would imply the AGN contribution to be 10 times larger in 
these two sources compared to the rest of the sample. 

Two other individual lines are worth mentioning; the [S IV] 
line, and the 18/um and 26/um [Fe H] Unes. The [S IV] hne 
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is present in most of the objects but is always weak. The 
line has an IP of 34.8 eV, below that of the [Ne III] line, 
and therefore should be easily excited. However, the photo- 
ionization cross-section of S^""" is lower than Ne^^, and one 
needs many photons to produce the [S IV] line. The IP of Fe 11 
is even lower, at 7.87 eV. It could therefore ionize within the 
photo-dissociation region (PDR). Lebouteiller et al. (2008) 
have shown that this ion correlates well with the [Ar II] line 
which is a good tracer of the ionized gas near the PDR. 

Fine structure hnes of the same species but from dif- 
ferent ionization stages can be used to trace the prop- 
erties (hardness and ionization parameter) of the ioniz- 
ing radiation field. The diamond symbols in Figure 6 
shows the ratio of the [S IV]10.5;(im (IP=34.8 eV) and 
[S III]18.7/xm (IP=23.3 eV) hnes plotted against the ra- 
tio of the [Ne III]15.5/im (IP=41.0 eV) and [Ne II]12.8Atm 
(IP=21.6 eV) lines for the starburst sample. As expected 
both ratios trace each other. Exceptions are Mrk266 and 
NGC 1365. These two outliers both have a strong AGN con- 
tribution to the neon lines, as measured by their >10 times 
elevated [Ne V]/[Ne II] ratio compared to the other sources. 
Excluding these two sources, the slope for the correlation be- 
tween [S IV]/[S III] and [Ne III/[Ne II] is 0.25 ±0.03. Our val- 
ues values fall within the range of values found by Verma et 
al. (2003) in their ISO study of starburst galaxies (including 
NGC 1365). Figure 6 also includes representative subsamples 
of Blue Compact Dwarf Galaxies (BCDs Hao et al. 2009), 
ULlRGs (Farrah et al. 2007), and galaxies from the Spitzer In- 
frared Nearby Galaxy Survey (SINGS Dale et al. 2009). The 
ratios in the SINGS galaxies studied by Dale et al. (2009) and 
in several star forming regions in irregular galaxies analyzed 
by Hunter & Kaufman (2007) follow the trend of the star- 
burst sample. The ratios in the starburst galaxies are lower 
than the BCD sample studied by Wu et al. (2008), and in fact 
Figure 6 shows only the lowest BCD ratios from their sam- 
ple. Hao et al. (2009) compare in detail our starburst sam- 
ple with similarly analyzed samples of ULIRGs, BCDs, and 
AGNs. They find that the starburst galaxies have the lowest 
excitation and are located in the lower branch of this diagram. 
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Fig. 9. — Abundance relations for those objects with a measured HI line in 
the IRS spectrum. These include 14 galaxies of which three have an AGN. 
Top panel: Neon and sulfur abundances. Bottom panel: Neon over sulfur 
ratio plotted against the neon abundance as an indicator of the metallicity. 
The grey box in the bottom panel indicates the range of solar values found 
in the literature (see §5.2). Note that the solar sulfur abundance ranges from 
1.4 to 2.1 X 10"^, which is significantly higher that the range of values in the 
starburst sample (top panel) and is therefore not plotted. 
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Fig. 10. — Neon abundance versus the hardness of the radiation field. 
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TABLE 5 

Elemental Abundances". 



Object S^+ZH S^+/H Ne+/H Ne++/H Fe+/H Fe++/H ICF S/H*" Ne/H Fe/H 
(10-«) (10-'') (10-*) (10-5) (10-«) (10-«) FelV (10-«) (10^) (IQ-*) 



NGC 1097 


3.22 


<0.22 


0.84 


0.68 




0.79 


0.00 


0.14 


0.11 


NGC 1222 


3.53 


2.09 


0.49 


2.58 




0.87 


0.37 


0.08 


0.43 


NGC 1365 


1.96 


1.47 


0.56 


1.15 




0.48 


0.25 


0.10 


0.19 


IC 342 


5.95 


0.15 


1.27 


0.37 




1.46 


0.03 


0.21 


0.06 


NGC 1614 


5.42 


0.78 


1.80 


2.19 




1.33 


0.14 


0.31 


0.36 


NGC 2146 


3.23 


0.18 


1.18 


0.82 




0.79 


0.03 


0.20 


0.14 


NGC 3256 


3.21 


0.17 


1.07 


0.64 




0.79 


0.03 


0.18 


0.11 


NGC 3310 


4.07 


1.53 


0.61 


2.98 




1.00 


0.27 


0.10 


0.50 


NGC 3556 


2.01 


<0.16 


0.39 


0.28 




0.50 


0.00 


0.07 


0.05 


NGC 4088 


1.53 


<0.14 


0.50 


0.16 




0.37 


0.00 


0.09 


0.03 


NGC 4194 


3.86 


1.08 


1.00 


1.56 




0.95 


0.18 


0.17 


0.26 


NGC 4676 


2.74 


0.23 


0.75 


0.60 




0.68 


0.04 


0.13 


0.10 


NGC 48 18 


4.87 


0.26 


1.39 


0.49 




1.20 


0.04 


0.23 


0.08 


NGC 7714 


4.32 


1.35 


0.60 


2.17 




1.06 


0.24 


0.10 


0.37 



Solai^ 



1.18 


0.56 


1.09 


3.54 


0.91 


1.90 


0.19 


0.16 




0.79 


0.15 


0.35 


0.40 


0.78 


1.58 


4.09 


0.75 


1.87 


0.07 


0.22 




0.91 


0.12 


0.29 


0.56 


0.31 


1.17 


2.30 


0.68 


1.02 


0.09 


0.08 




0.50 


0.12 


0.17 


1.47 


1.34 


1.07 


6.56 


1.31 


3.01 


0.24 


0.38 




1.46 


0.22 


0.62 


2.23 


2.94 


1.20 


6.04 


2.02 


6.21 


0.36 


0.82 




1.34 


0.35 


1.18 


0.64 


0.37 


1.09 


3.57 


1.26 


1.10 


0.10 


0.10 




0.79 


0.21 


0.21 


0.71 


0.87 


1.13 


3.55 


1.13 


1.78 


0.12 


0.25 




0.79 


0.19 


0.37 


1.08 


0.90 


1.38 


4.63 


0.91 


2.73 


0.18 


0.25 




1.03 


0.15 


0.43 


0.21 


0.19 


1.13 


2.21 


0.42 


0.45 


0.04 


0.06 




0.50 


0.07 


0.09 


0.30 


0.18 


1.06 


1.68 


0.52 


0.51 


0.05 


0.05 




0.37 


0.09 


0.10 


0.69 


1.06 


1.25 


4.35 


1.16 


2.19 


0.11 


0.30 




0.97 


0.20 


0.41 


0.65 


0.00 


1.00 


3.04 


0.81 


0.65 


0.11 


0.00 




0.68 


0.14 


0.11 


2.09 


1.99 


1.08 


5.38 


1.44 


4.43 


0.35 


0.58 




1.21 


0.24 


0.92 


0.71 


1.22 


1.43 


4.89 


0.82 


2.76 


0.12 


0.34 




1.08 


0.14 


0.46 








1.4-2.1 


0.7-1.2 


32.4 



"Abundances are given in number with the uncertainties indicated in the row below each measurement. 
"^Includes 10% of S III to account for the contribution of S II (see §5.2). 

"^Range of solar abundances taken from Grevesse & Sauval (1998); Feldman & Widing (2003); Asplund et al. (2005, 2006). 
See §5.2 for explanation. 



Compared to the Starburst and ULIRG samples of Verma et 
al. (2003) and Farrah et al. (2007) the excitation level in our 
sample is on the low and middle end, with the highest ratios 
being [S IV]/[S III]-0.4 and [Ne III]/|Ne II]--1, and the Verma 
et al. (2003) and Farrah et al. (2007) samples reaching val- 
ues as high as 10 for both ratios. This lower ionization in our 
sample may indicate that there are fewer massive stars or that 
they are older Interestingly, the three sources with the high- 
est [Ne III]/[Ne II] ratios in our sample (>0.6) are not those 
with an AGN component, but pure starbursts. NGC 1222, 
NGC 3310 and NGC 7714 all are found on the 'BCD branch' 
of the [Ne V]/[Ne H] versus [Ne ni]/[Ne H] diagnostic diagram 
ofHaoetal. (2009). 

Groves et al. (2008) smdied the relation between the well 
known mid-infrared ionization diagnostic with the ground ac- 
cessible ratio of [S IV]/[Ne II]. They found a good correla- 
tion between the ratios in their sample, which included star- 
burst galaxies, BCDs, ULIRGs, AGNs, Planetary Nebulae, 
and Galactic and extra-galactic H II regions. We confirm this 
relation (Figure 7), and find a slope value of 0.68 ±0.08 (ex- 
cluding NGC 4945), in good agreement with the slope found 
by Groves et al. (2008) (0.65 for the starburst sub-sample and 
0.81 for their fuU sample). NGC 4945 has the strongest sili- 
cate absorption feature of the sample (Brandl et al. 2006) and 
this attenuates the [S IV] line flux relative to the [Ne H] line 
which is exactly what is seen in Figure 7. 

Ho & Keto (2007) illustrated how the neon lines in the mid- 
infrared spectra can be a powerful indicator to derive the star 



formation rate in star-forming galaxies. They showed that the 
sum of the luminosities of both lines correlate well with the 
infrared luminosity in over 5 orders of magnitude in luminos- 
ity, and proposed a calibration between the luminosity of the 
lines and the star formation rate. Figure 8 shows the relation 
between the neon fines and infrared luminosity. Our data cov- 
ers a narrow range of infrared luminosities, but it can be seen 
that a correlation is present. The slope in the figure has a value 
of 1 . 1 8±0.20, agreeing within the errors with that reported by 
Ho & Keto (2007) (0.98±0.069). 

5.2. Abundances 

Infrared lines offer many advantages over optical or ultravi- 
olet (UV) fines for measuring elemental abundances, the most 
important of which is their small dependence on the electron 
temperature (e.g. Bernard-Salas et al. 2001; Wu et al. 2008). 
This reduces the uncertainty in the abundance estimates when 
there are electron temperature variations, which is the case in 
integrated spectra of galaxies, or when the electron tempera- 
ture cannot be measured. IR lines are also less affected by ex- 
tinction. Optical and UV studies also rely on the often rather 
uncertain, ionization correction factors (ICFs) to account for 
unobserved stages of ionization. The direct measure of im- 
portant stages of ionization of some elements in the infrared 
greatly reduces, or completely avoids, the use of ICFs. 

The HI (6-7) line at 12.37/im is detected in 14 of the 
stabursts. This allows us to measure the abundances of neon, 
sulfur, and (tentatively) iron relative to hydrogen. The abun- 
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dances are derived adopting an electron temperature (T^,) of 
10000 K and a gas density of 100 cm"-'. These are values 
typical of H II regions. We derived abundances for a range of 
Te (7500-12 500 K), and used the differences to obtain the un- 
certainties in the measurements. These uncertainties amount 
to around 20% in the total abundances. Using these parame- 
ters and solving the equation of statistical equilibrium for the 
levels of the atoms, the ionic abundance can be derived (see 
eq. 1 Bernard-Salas et al. 2001). These values are reported in 
Table 5. Note that we cannot derive the density from the pair 
of [S III] (18.7^m/33.4/im) or [Ne III] (15.5/im/36.0^m) Hues 
because the dependence of these ratios on the density starts 
above 3000 cm"^ for the sulfur ratio, and higher for the neon 
ratio. Compared to our reported values, an electron density 
of 1000 cm"^ (on the high side for H 11 regions) will produce 
abundances that are no more than 1.7% higher for Neon, 13% 
higher for Sulfur, and about 80% lower for Iron. The latter is 
due to the electronic configuration of iron. 

In starburst galaxies the most important stages of ioniza- 
tion of neon for abundance determination are Ne"^ and Ne^ 
which are measured in the IR. Thus, no ICFs are needed. 
Among the 14 objects with a HI detection, Ne"*""" is present in 
NGC 1365 and NGC 4194 (Table 5), but this likely originates 
from the AGN given the hard radiation field needed to break 
the high ionization potential (IP) barrier of 97.1 eV to pro- 
duce this ion. This applies as well to the Ne-'"^ (IP of 63.5 eV, 
also above the double ionized helium cutoff). In addition, the 
contribution of Ne"*"^ to the total neon abundance is less than 
1%. The [Ne V]14.3Aim/[Ne II]12.81/im ratio in NGC 1365 
and NGC4194 is 0.13 and 0.018 respectively. These ratios 
correspond in the Farrah et al. (2007); Spoon et al. (2009, in 
prep) diagrams to an AGN contribution of ^10% and <1% 
respectively. This implies that the abundances we derive for 
NGC 4194 is indeed representative of the starburst while in 
NGC 1365 our abundances may be ~10% too high. The to- 
tal sulfur abundance is derived adding the contribution from 
S^^ and S^^. No S^ is measured in the infrared but its contri- 
bution to the total sulfur abundance is usually small. Martin- 
Hernandez et al. (2002) estimated a contribution of about 15% 
of S"^ in the Galactic ultra-compact HII regions they studied. 
Vermeij & van der Hulst (2002) analyzed a sample of Magel- 
lanic Cloud H II regions using IR data, and combined optical 
measurements of [S II] lines. They derived an ionic abun- 
dance of S^ between 2 to 20% of the total sulfur abundance, 
with ~8% being the average. Using the photo-ionization code 
of Stasinska (1990) we have created a grid of models at typi- 
cal conditions of H II regions; stellar temperatures considered 
between 37 500 to 55 000 K, and electron density of 10 cm""*, 
and a filling factor ranging from 0.2 to 1.0. These models pre- 
dict the contribution of S^ to be between 7 and 15% of that 
of S"^ (which is the dominant contributor to the total sulfur 
abundance). We have thus adopted a 10% contribution of S"*" 
to the sulfur abundance (ICF=1 . 1) to account for the contribu- 
tion of this ion in our sample. While we have measured Fe^ 
and Fe"^ from our Spitzer data, Fe^^ may be an important con- 
tributor (up to half) to the total iron. Its contribution has been 
estimated using the photo-ionization grid models by Stasinska 
(1990), and taking advantage of the similarity in IP between 
the [Ne III]/[Ne II] and [Fe IV]/[Fe III] line ratios to scale the 
conttibution of this ion in each object. At can be seen from 
Table 5 the ICF varies from a few percent to a factor 2. Given 
the importance of this ion to the total iron budget in those ob- 
jects with a high ICF, the iron abundances we report should 
be interpreted with care. The uncertainty for this element in 



Table 5 only takes into account the uncertainty in the Fe^ and 
Fe"^ ionic abundance. 

In Figure 9 (top), the neon abundance of our starburst sam- 
ple correlates reasonably well with the sulfur abundance. This 
is expected because they are both a-elements and are be- 
ing produced in the same type of stars. This finding was 
also pointed out by Wu et al. (2008) in their BCD sample 
and it contradicts earlier findings Verma et al. (2003) who 
based on ISO data did not find such a correlation between the 
abundance of these two elements for the starbursts they stud- 
ied. In our sample the neon abundances range from slightly 
supersolar (IC342, NGC 1614, NGC 48 18) to several times 
lower than the solar value (see Table 5). Note that the so- 
lar abundances of certain elements have been subject to many 
changes during the last decade (Pottasch & Bernard-Salas 
2006; Bemard-Salas et al. 2008; Lebouteiller et al. 2008), 
and thus, we resort to compare to the range of values found 
in the literature (0.69-1.2x10"^). There is however increas- 
ing evidence that the higher solar neon abundance reported 
agrees better with the findings of studies in Galactic HII re- 
gions and planetary nebulae. The sulfur abundance is also 
lower than solar (more so than neon, with solar sulfur val- 
ues ranging from 1.4-2.1 x 10"^), which is in agreement with 
many previous studies. Available optical-derived abundances 
of oxygen by Pilyugin et al. (2004); Pastoriza et al. (1993); 
Gonzalez-Delgado et al. (1995) in IC342, NGC 3310, and 
NGC 7714 (7.1, 1.5,1.8x10"^ respectively) confirm the trend 
of our IR-derived neon abundances, with IC 342 supersolar, 
and the other two galaxies having a 2-3 times lower abun- 
dance than the Sun (adopting 4.9 x 10"^ for the solar oxygen 
abundance). 

The bottom panel in Figure 9 displays the Ne/S ratio versus 
the neon abundance. The range of solar values is indicated by 
the grey box in the figure. The Ne/S ratio varies from 18 to 
40. This range of values is typical for other infrared derived 
abundances in starburst galaxies, HII regions and even plan- 
etary nebulae (Rubin et al. 2007, 2008; Bernard-Salas et al. 
2008). It is overall higher than the ratio found by Wu et al. 
(2008) in their abundance analysis of BCDs, where they find 
values from 10-20. The Ne/S ratios are higher than the solar 
Ne/S ratio, revealing a discrepancy between the neon and sul- 
fur solar values, and which has now been widely reported (e.g. 
Marigo et al. 2003; Henry & Kwitter 2004; Rubin et al. 2007; 
Bernard-Salas et al. 2008; Wang & Liu 2008, and references 
therein). In Figure 10 the neon abundance is plotted versus 
the hardness of the radiation field. The sample is small but 
seems to indicate that a hard radiation field is found only for 
sources with a low metallicity which ties well with the com- 
mon knowledge that at low stellar metallicity, line blanketing 
is diminished resulting in harder stellar radiation fields. 

5.3. PAHs 

In Figure 1 1 the PAH luminosity for the different bands is 
plotted. We chose to compare the PAH features in the spectra 
to the integrated emission of the 1 1.0+1 1.3/xm complex be- 
cause it is the strongest in this wavelength region (panels 'a' 
through 'j'). Panels e, h, and j show a very good correlation 
of the 1 1.3/im complex with the relatively strong PAH bands 
at 12.7, 16.4, and the 17.0/im PAH complex. The 17.0yum 
complex was discovered by Van Kerckhoven et al. (2000) as a 
weak feature in 750 spectra and attributed to PAHs by labora- 
tory analysis (Leger & Puget 1984; Allamandola et al. 1985). 
Smith et al. (2007) corroborate this using spectra from the 
SINGS galaxies and showing good correlation between the 
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Fig. 11. — Correlation between the PAH bands. 
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TABLE 6 
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Fig. 12. — Top panel: Sum of PAH luminosities versus the H2 emission. 
Bottom panel: the 11.3/im PAH equivalent width with respect to the hardness 
of the radiation field as given by the [Ne III] 15.5//m/[Ne II] 12.8/um ratio. 



Strength of the and the 1 1.3/im PAH complexes. Our 

sample supports this finding (panel 'j'). Other bands also cor- 
relate well (e.g 11.0, 12.0, 16.7/im). The correlation of the 
weaker bands (10.6, 13.5, 14.2/xm) is good but shows more 
scatter. The 16.7/zm feature is listed by Smith et al. (2007) as 
part of the 17.0/^m complex but in our higher resolution spec- 
tra we can separate it from this complex. Smith et al. (2007) 
also report a PAH feature at 15.9//m. This feature is not 
present in our high S/N template spectrum but may be present 
in NGC 3079, NGC 3628 and NGC 4676, albeit weakly. The 
unidentified feature at 10.75//m (panel 'b') shows a good cor- 
relation with the 1 1.3/im complex. In the bottom panels (T 
and 'm') we plot this feature with respect to other weak PAHs. 
Although there is some scatter, they do seem to correlate. It 
is therefore likely that the 10.75/um feature is either a PAH or 
at least its carrier is carbonaceous in nature. Panel 'n' shows 
a very good correlation between the satellite features in the 
1 1.0 and 12. 0/im bands, and the last panel displays the equally 
good correlation between the lljim PAHs. In summary, the 
strengths of the PAH features scale with each other. 

Figure 12 displays the relation between the PAH features 
and the emission lines. In the top panel, the sum of the PAH 
luminosities is plotted against the sum of the H2 lines. Both 
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^Based on the 0-0 S(l) and S{2) transitions. 



are produced in the PDR, albeit in different regions, and their 
fluxes correlate well. In the bottom panel the equivalent width 
(EW) of the 1 1 .3 /xm PAH'^ is plotted against the [Ne ni] over 
[Ne II] line ratio. There is no correlation, indicating that there 
is no dependence of the PAH with the hardness of the radia- 
tion field. In fact, the EW of the 1 1.3/im shows little variation 
in Figure 12. In AGNs, the PAH EW is seen to decrease with 
harder radiation field as the PAH are further ionized or de- 
stroyed, but the PAH EW is flat for soft radiation fields as 
those in Starbursts (Lebouteiller et al., in prep.). 

Another characteristic of the PAHs is their profile, which 
has been found to change between and within sources. This 
led Peelers et al. (2002) to establish a classification based on 
the peak position of the different bands. Figure 13 shows the 
continuum subtracted spectra for all the objects in the sample 
in the regions of interest (where the emission lines have been 
removed). The grey area in the figure represents the template 
spectrum. As opposed to the template shown in Figures 3 
and 4, this template has been normalized to the sum of the 
PAHs (and not the sum of the PAHs plus lines), but the are 
both very similar. This PAH continuum subtracted spectrum 
of the template starburst is also available in electronic format 
(tab26.txt). The fact that the spectrum of each object (thin 
black line in the figure) follows that of the template (grey area) 
indicates there is hardly any variations in the profiles. 

The relative strengths of the features is also the same 
for most of the object although there are a few variations. 
The most striking is 12.7/im PAH in NGC 4945 which is 
stronger than the 11.3/xm PAH. The reason for this is the 
stronger silicate dust extinction on the 11.3/um PAH than the 
12.7/imPAH. This band is also strong (relative to the 1 1.3/im 
PAH) in those starburst with AGN signature (e.g. NGC 660, 

For simplicity we use the equivalent width of the combined 11.0 and 
11.3 /im components. This is valid as both bands correlate well with each 
other (Figure 1 1). 
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Fig. 13. — Continuum subtracted spectra of the starburst sample. The thin black line represents the spectrum of the corresponding object and the grey shaded 
area that of the template (with a white edge to distinguish the grey areas in the overlapping 11.3//m band). The fine-structure line emission has been removed to 
better show the PAH profile in the starburst spectra. The spectra have been normaUzed to the total PAH emission and ordered according to their 1 1.3 /im PAH 
strength (strongest at bottom), and the 15-18;im flux density has been multiphed by a factor 5 to better display the PAHs in that region. 



NGC 3079, NGC 3628) but also in a few objects with no AGN 
(e.g. NGC 253). The 17.0/im complex does not always follow 
that trend. Despite these differences the profiles are near akin. 
As we saw from previous sections the metallicities are similar 
(within a few factor of each other), and the radiation fields are 
low in most of them. Therefore, it is possible that the invari- 
abiUty in the PAH spectra of the starburts in our sample are 



probably due to the similar conditions in these galaxies. 

5.4. Molecular Hydrogen 

We detect three pure-rotational transitions 0-0 S(2), S(l), 
and S(0) of H2 at 12.28, 17.03, and 28.21 fim respectively. 
The S(l) transition at 17.03 fim is the strongest of the three. 
Figure 14 shows that, while there is some scatter, there is a 
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Fig. 14. — Luminosities of the [Si II] line versus the sum of the three H2 
rotational lines measured in the IRS spectrum. For this plot the fluxes of the 
H2 S(2) and S(l) lines have been scaled based on the SH to LH continuum 
(see text for further details). 



good correlation of the H2 lines with the [Si II] line. This is 
expected since both originate in the PDR. In fact, |Si II] is 
an important fine-structure cooling line of the PDR. We note 
that in order to combine the three lines we used a scaling fac- 
tor between the SH and LH modules determined by matching 
the continuum in the overlap region in both modules around 
19.5 /xm. This is not ideal since no off observations were 
subtracted and it contributes to the scatter present in the fig- 
ure. The ratio of the sum of the H2 lines over the [Si II] line 
is 0.22, with values ranging from 0.04 to 0.5. This is very 
similar to the values found by Roussel et al. (2007) in their 
starburst (SINGS) sample. 

The molecular hydrogen lines can be used to give informa- 
tion on the conditions of the warm component in the PDR. 
There is a relation between the density and temperature as- 
suming local thermal equilibrium and using the Boltzman 
equation (van den Ancker et al. 2000; Bemard-Salas & Tie- 
lens 2005). Using this relation and the lines in the SH module 
(to avoid scaling effects), we have derived the excitation tem- 
perature (Tex), column density N(H2), and warm molecular 
mass for the sample. These are given in Table 6. The aver- 
age excitation temperature and column density are 360±64 K 
and 2x 10^" cm"^ respectively, and the warm molecular mass 
ranges from several 10"* to SxIO^Mq, with an average mass 
of 1 . 1 X 1 0''M0 . We stress that these numbers correspond only 
to the region covered by the SH sht and are based on the fluxes 
for just the S(l) and S(2) Unes. 

Recently, Roussel et al. (2007) studied the warm molecu- 
lar hydrogen emission in the nuclear region of a sample of 
SINGS galaxies. The range in areas covered in their observa- 
tions (60 pc to 3.8 kpc) is very similar to the range covered 
in our study (see §2). Their column densities compare well to 
those found for our starburst galaxies, but their excitation tem- 
peratures are on average cooler by ~ 150 K, and their masses 
significantly larger (10^-10^ Mq). This is expected because 
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the SINGS galaxies have a much lower star formation rate per 
unit volume compared to our starbursts, and since they are 
less luminous, a larger fraction of their molecular Hydrogen 
is colder. 'Cold' molecular hydrogen is traced by the S(0) line 
which is not used in our analysis as it falls in the LH sUt, but it 
is used in the SINGS sample where aperture size differences 
do not exist. Hence, our computed H2 masses are lower than 
those found for the SINGS sample, but this is hkely an artifact 
of the absence of the S(0) line in our calculations, rather than 
a real difference. 

We find no relation between the excitation temperature and 
sources with and without AGNs. Higdon et al. (2006) an- 
alyzed Spitzer data of a sample of ULlRGs. Their average 
excitation temperature (336±15 K) agrees with that of our 
starburst sample, but again their warm molecular masses are 
several hundred's times larger (2x10**Mq). It is important 
to note that most ULlRGs fit in the SH aperture while in the 
starbursts we are only measuring the central region. Still, this 
difference may point as well to a larger fraction of mass being 
colder in the Starburst sample. 

6. CONCLUSIONS 

1. We have presented high resolution spectra (R~600) of 
a sample (24) of nearby starburst galaxies. The spec- 
tra are dominated by emission from forbidden fine- 
structure lines and polycyclic aromatic hydrocarbons 
features. 

2. [Ne V] is detected in six of the eight starbursts with a 
previously known AGN component. This line is also 
seen in NGC 4194. 

3. Several PAH features have been measured in the spec- 
tra, some of which are weak and are rarely seen in the 
spectra of starburst galaxies and HII regions. 

4. An emission feature around 10.7/im and an absorption 
feature at 13.7/xm are detected for the first time. Their 
identification is unknown, but the strength of the feature 
at 10.7/im correlates with that other PAH features. 

5. The average conditions of the warm gas as traced by 
the S(l) and S(2) pure rotational molecular hydrogen 
lines indicates an average excitation temperature of 
355±60 K and a density of 2x 10^° cm'^. The warm 
molecular mass ranges from several 10* to 4 x 10* Mq . 

6. We provide the combined average spectrum of the 15 
'pure' starburst galaxies of our sample, which is likely 
the highest signal to noise mid-infrared spectrum of a 

10" L0 starburst galaxy available from Spitzer. This 
template is representative of a large range in luminos- 
ity and IR color, and provides a good indication of the 
spectral properties of starburst galaxies at high redshift. 
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